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I ' Abstract. We introduce an approximation to calculate the gravitational radiation pro- 
^ ■ duced by the collision of true-vacuum bubbles that is simple enough to allow the sim- 
^ ■ ulation of a phase transition by the collision of hundreds of bubbles. This "envelope 
approximation" neglects the complicated "overlap" regions of colliding bubbles and fol- 
lows only the evolution of the bubble walls. The approximation accurately reproduces 
j_| ■ previous results for the gravitational radiation from the collision of two scalar-field vac- 
^ uum bubbles. Using a bubble nucleation rate given by F = Foe^*, we simulate a phase 

transition by colliding 20 to 200 bubbles; the fraction of vacuum energy released into grav- 
ity waves is -E'cwZ-E'vac = 0.06(if//9)^ and the peak of the spectrum occurs at Wmax = 1-6/9 
{H"^ — 8ttGp/3 is the Hubble constant associated with the false- vacuum phase). The spec- 
trum is very similar to that in the two-bubble case, except that the efficiency of gravity- 
wave generation is about five times higher, presumably due to the fact that a given bubble 
collides with many others. Finally, we consider two further "statistical" approximations, 
where the gravitational radiation is computed as an incoherent sum over individual bub- 
bles weighted by the distribution of bubble sizes. These approximations provide reasonable 
estimates of the gravitational-wave spectrum with far less computation. 
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I. Introduction 



The cosmic background of gravitational radiation provides a unique probe of the early- 
universe. Unlike electromagnetic radiation, gravity waves propogate virtually unimpeded 
since the Planck epoch, providing an unmodified record of cosmic events. Possible cosmo- 
logical sources include the thermal background (the graviton analogue of the microwave 
background), inflation [1], cosmic strings [2], a pregalactic star population [3], and phase 
transitions [4]. In particular, strongly flrst-order phase transitions are among the most 
promising of all these sources [5] : The energy released in gravitational waves can approach 
1% of that in the ambient thermal bath. 

In a recent paper we initiated a detailed investigation of gravity wave production 
from strongly first-order phase transitions by calculating the radiation from two colliding 
vacuum bubbles [6]. Beginning with a scalar field configuration corresponding to two 
bubbles nucleated simultaneously and far apart, we used the Klein-Gordon equation to 
evolve the scalar field for a time r comparable to the initial bubble separation. (In realistic 
phase transitions the duration of the transition is comparable to the typical separation of 
nucleation sites [7].) From the scalar-field conflguration, we calculated the stress-energy 
tensor and, in the linearized gravity approximation, the energy spectrum of radiated gravity 
waves. The pair of vacuum bubbles radiates efficiently: the fraction of energy that goes 
into gravity waves is 



where E'vac is the total energy liberated by the two vacuum bubbles, r is the total time 
of the bubble evolution, expected to be of order 0.01 to 1 of and H'^ — 87rGpvac/3 

is the Hubble parameter associated with the vacuum energy. The spectrum of radiation 
peaks at a frequency 



These results imply that vacuum bubble collisions can indeed be potent sources of 
gravitational radiation. However, our work has its limitations, the foremost being the use 
of a time cutoff in the simulation to model the end of the phase transition. Specifically, we 
smoothly ramped the scalar field gradients (the source of gravitational radiation) to zero 
after a time r. While such an ad hoc prescription greatly simplifies the problem and is 
probably reasonable, it clearly neglects multi-bubble effects. The motivation for colliding 
many bubbles is to model more realistically a phase transition. 

A direct attack on the many-bubble problem employing scalar-field evolution is numer- 
ically infeasible. The two- bubble problem was made tractable by exploiting the 0(2, 1) 
symmetry possessed by the space-time of two vacuum bubbles, which makes the scalar field 
evolution effectively one, rather than two, dimensional. Even the slight generalization to 
a pair of bubbles nucleated at different times, a situation still having rotational symmetry 
about the axis connecting the two bubbles, proved nearly impossible. The most general 
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case of many bubbles in three dimensions has no symmetries and is beyond present com- 
puting capabiUties. The problem is difficult numerically because of two disparate scales. 
The bubble wall thickness at nucleation is small compared to the size of the bubble at 
collision; moreover, the bubble wall becomes thinner due to Lorentz contraction as the 
bubble expands. 

To proceed further requires dispensing with the detailed dynamics of the scalar field. 
The results of our two-bubble simulations suggested an elegant approximation. The spec- 
trum and amount of gravitational radiation depended only on the gross features of the 
bubble collisions: the vacuum energy and the size of the bubbles at the end of the phase 
transition {i.e., the cutoff time). Even though the field dynamics after a bubble collision 
are quite intricate, the overall contribution to the radiation from the complicated small- 
scale motions adds incoherently and is subdominant. This prompted us to consider an 
"envelope approximation": the bubbles are treated as infinitely thin, and in the regions 
where bubbles overlap, the bubble wall is completely ignored. Only the envelope of the 
evolving bubble network is considered. As we shall discuss, the envelope approximation 
very accurately reproduces our previous results for two colliding bubbles and allows us 
to model a phase transition with the collision of hundreds of bubbles. When applied to 
a phase transition where the bubble nucleation rate increases exponentially with time, 
F oc exp pt, the fraction of vacuum energy liberated in gravitational waves is found to be 



or about five times the efficiency estimated from the collision of two bubbles (for such a 
nucleation rate, the duration of the transition is r ~ few P) . 

The paper is organized as follows. The next Section of this paper describes the envelope 
approximation, with detailed comparisons to the previous two-bubble results obtained from 
scalar- field evolution [6] . In Section III we first review some pertinent aspects of bubble 
nucleation theory [7] , and then present our numerical results for the gravitational radiation 
from large numbers of colliding bubbles. In Section IV, we present another approximation 
which treats the production of gravitational waves in a statistical sense: as the incoherent 
sum of radiation from individual bubbles, weighted by the distribution of bubble sizes. 
Many of the results in Section III can be reproduced by this simple recipe, and with far 
less computation. We finish with a summary of our work and some concluding remarks. A 
review of tensor spherical harmonics and some auxiliary formulas used in Section IV are 
relegated to an Appendix. A summary of the results of this paper and our earlier work 
[6], as well as the application of our results to cosmological phase transitions, is presented 
elsewhere [8]. 
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II. Envelope Approximation 

(a) Review of vacuum bubbles 

We consider a real scalar field (fi with a potential possessing two non-degenerate local 
minima: 

jC=-d^^^d,^-Vicp). (4) 

Throughout we use a metric with signature (H ). The exact form for the potential is 

not important, but where needed we use 



The dimensionless number e measures the degree of symmetry breaking between the two 
minima near ±(/?o- The relative minimum corresponding to +ipo is the "false vacuum," 
while the global minimum corresponding to —<f>o is the "true vacuum." The vacuum energy 
density is defined as the difference in energy density between the true and false vacua; 
here, pvac — 2eA<^Q. The height of the potential barrier between the two vacuum states is 
~ X(Pq/S. The relevant features of the potential are that it possesses two inequivalent local 
minima differing in vacuum energy by Pvac? and that the height of the barrier between the 
two minima is large enough so that the false vacuum decays via quantum tunnelling. 

Classically, the false-vacuum state is stable, but quantum effects cause its decay to 
the true-vacuum state. This decay proceeds via the quantum nucleation and expansion of 
bubbles of the true- vacuum phase which spontaneously appear from the false- vacuum state. 
Coleman has shown that the bubble with minimum action is 0(4)-invariant in Euclidean 
space [9]; the initial bubble profile is obtained by analytically continuing to Minkowski 
space and taking the t = time slice. The vacuum bubble then evolves according to the 
Klein-Gordon equation and has 0(3, 1) symmetry; i.e., the scalar field <^ is a function only 
of the quantity — x^ —y^ — z^. The energy difference between the true and false vacuum 
phases creates an effective outward pressure on the bubble wall, causing it to expand with 
constant acceleration. For our purposes, the important aspects of bubble dynamics are 
that the expansion speed rapidly approaches the speed of light, and that the false- vacuum 
energy liberated becomes kinetic and gradient energy of the bubble wall [10]. 

The bubble's symmetry allows us to quantify the above statements while deriving a 
result which will be useful in the envelope approximation. First, the 0(3, 1) symmetry 
immediately implies that the position of the bubble wall is given by 



where Rq is the initial radius of the bubble and x^aii denotes a fiducial point within the 
bubble wall. Next, consider the stress-energy tensor associated with the expanding bubble: 
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The energy density in the scalar field is given by the time-time component of the stress 
tensor: 
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where we have used the spherical symmetry of the bubble solution. The bubble's 0(3, 1) 
symmetry can be used to write the the energy of the bubble wall at any time t as 



E{t) pa 471 J drr 
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where s = ^/r"^ — t'^ and f{s) = ip{r, t = 0) is the profile of the initial bubble solution. Note 
that for < t^, d(p/ds ~ so the integral is zero. The two terms represent the kinetic 
and gradient energy of the bubble wall. We neglect the potential energy term inside the 
bubble wall as it rapidly becomes unimportant as the bubble wall gets thinner; we have 
not included the false-vacuum energy outside the bubble, as we are only interested in the 
energy liberated by the bubble. In Fig. 1, the kinetic and gradient energy for a bubble are 
shown as a function of time. After a short time, the energies become equal; each is half of 
47rpvac^"^/3, the total vacuum energy liberated by the bubble. From Eq. (9b) it is simple 
to see that gradient and kinetic energies are equal and increase as for t ^ Rq: dip/ds is 
only nonzero when s is close to Rq] when t ^ Rq each term in the integrand approaches 
{d(f/ds)H^/s. 

(b) Envelope approximation 

As in our previous paper [6], we compute gravity- wave production in the linearized 
gravity approximation, valid for bubble sizes less than (recall = 87rGpvac/3). The 
energy radiated in gravitational waves can be expressed in terms of the Fourier transform 
of the spatial components of the scalar-field stress-energy tensor. Further, in computing 
Ty(k, a;) we may neglect the Cgij piece as it is a pure trace and does not act as a source 
for gravitational radiation. Thus the fundamental quantity is 



T 



1 /"°° f 
(k, co) = —J dt e*'^* J d^xdiifdjipe- 



(10) 



where k is a unit wave-vector. As before, we adopt Weinberg's unusual normalization 
convention for the Fourier transform [11] . The scalar-field configuration of interest is that of 
many colliding vacuum bubbles. In the envelope approximation we assume that the overlap 
regions where bubbles have expanded into one another do not contribute substantially to 
the gravitational radiation, and exclude these regions from the spatial integration (see 
Fig. 2). We can then break up the integral into integration regions, one surrounding each 
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nucleation site and extending out to the bubble radius. Equation (10) becomes 



Tij(k,uj) ^ — 



l>oo 


r ^ r 






Jo 





dn'e 



I ia;k-x 



(11) 



where is the number of bubbles, Sn is the portion of the surface of bubble n that remains 
uncoUided at time t, and the primed spherical coordinates are chosen independently around 
the center of each bubble. We have also assumed that the wall thickness is small compared 



to uj ^: u}/S.r <^ 1. This means e" 



is essentially constant across the bubble wall and 



can be factored out of the r-integral. In practice, for the frequencies of interest, this is 
always an excellent approximation. 

Next we use the fact that each bubble is spherically symmetric around its center, 
so that is independent of the angular variables. Dropping the primes for notational 
convenience, the stress tensor components become 
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where x^ is the nucleation site of the nth bubble and Xi is the iih component of a unit 
vector pointing from x„ in the direction dO,: 

X = sm.9 cos 0, y = sin 6 sin 0, z = cos 9. 

Above we showed that the kinetic and gradient energy associated with a bubble wall are 
equal after a small amount of bubble expansion. Using Eq. (9a), for each bubble we have 

2 



J drr^ ^ 
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where the bubble radius R{t) ~ t. This is a very good approximation by the time a bubble 
has doubled in size; cf. Fig. 1. Substituting for the radial integral in Eq. (12) 

leads to 



dQ e 



where t„ is the nucleation time of bubble n. 
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The total energy radiated in gravity waves is given in terms of Tjj(k,a;) by [11] 



dE 
dujdVL 



= 2Gu;^Aij,irnik)T*Jk,u;)Tirn{k,u;) 



(15) 



where A^ij^im is the projection tensor for gravity waves, 
■^ij,Zm(k) = SiiSjjji — 2kjk^5j; + ikjkjk/k-^ 
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(c) Scaling properties 

From Eqs. (14) and (15), two important scaling relations are evident. First, the total 
radiated energy is explicitly proportional to Pvac J^^^ ^ '^^ found previously in the two- 
bubble case [6] . Second, since the bubbles expand at essentially the speed of light and to a 
good approximation have zero thickness and zero initial size, the problem has no intrinsic 
length/time scale. Making the transformation t — > ^t, x — > 7X, we find the following 
scaling properties: 

u) ^ —; (17a) 
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The length/time scale is set by the average separation between bubble nucleation sites. 

Equations (17) show that the total energy radiated from a volume containing a fixed 
number of bubbles will vary with the fifth power of the mean bubble separation. As 
we will show in the next Section, these scalings can also be expressed in terms of the 
bubble nucleation rate, since the typical separation of nucleation sites is determined by 
the nucleation rate. For the two-bubble case, the above scalings were found previously to 
hold to very good accuracy [6] , where the relevant scale is the total evolution time. These 
scalings have great practical importance as they allow us to apply the results of a single 
numerical simulation to any phase transition with a bubble nucleation rate of the same 
functional form. In particular, we shall use the functional form F oc e^* — a nucleation rate 
that increases exponentially with time — in which case the length/time scale is just 

(d) Two bubbles, quadrupole approximation 

In order to determine the accuracy of the envelope approximation, we compare the 
results it gives for two bubbles to the results calculated previously with the exact scalar field 
evolution [6]. As a warm-up, we begin with the conventional "quadrupole approximation," 
corresponding to the limit k • x ^ 0: 



T,,(k,a;) = ^ / dte'^' J2(t-t^f / 
Stt Jo Js 



(18) 



Note that because our source is not small compared to the wavelength of the radiation, this 
approximation does not correspond to the quadrupole term in the multipole expansion; 
see Ref. [6] for detailed elucidation of this point. The quadrupole approximation simplifies 
the requires calculations, and in the limit a; — > gives the correct result; thus we use it as 
a starting point for our comparisons. 

Consider two bubbles of negligible initial size nucleated simultaneously at t = on the 
^-axis at ±d/2. The bubbles will first "kiss" at t = d/2. Define cos a = d/2t for t > d/2 
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and ct = for t < (i/2; ct is the angle excluded from the angular integration because of 
bubble overlap in the envelope approximation. Then 
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Using the spherical coordinates defined in Eq. (11), 



Id/2 

and 
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T,, = ^ r dte'^' (t' + ^] Cit), (206) 



9 Jd/2 

where we have included a time cutoff function C{t) which decreases smoothly from 1 to 
on a time scale r ~ 0{d). As mentioned earlier, the cutoff was introduced in our previous 
two-bubble calculations to model the completion of the phase transition; C{t) is discussed 
in detail in Ref. [6]. 

Rotational symmetry around the 2;-axis implies the off-diagonal components of the 
stress tensor are zero, and must be of the form 

T,j = D{u;)Sij + A{u;)5i,5j,. (21) 

The first term, being a pure trace, does not contribute to gravitational radiation. The 
second term is 

'd^^ d 
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dte'^' — - -t^ C{t). (22) 



Substitution into Eq. (15) gives 



G'a;^|A(a;)|^sin^^; ^ = ^Gu^\A{uj)\\ (23) 



dwdQ duj 15 

The comparison of the envelope approximation with the previous calculations using 
scalar-field evolution is shown in Fig. 3. (For reference, in the scalar-field evolution case 
we used r/d = 1.2 and a gaussian roll-off in the final 10% of the evolution time for 
C(t); see [6].) The features of the spectrum are reproduced remarkably well. The overall 
normalization of the envelope approximation is high by 20 percent. It may seem strange 
that the total power radiated is higher from an approximation which neglects a chunk of the 
source. However, when two bubbles collide, a refiected wave begins to propogate outward 
from the point of collision. This wave takes the approximate shape that the colliding 
portion of the bubbles would have had they not collided. By neglecting the interaction 
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region, we actually make a given bubble less spherical, and hence increase the amount of 
radiation. 

(e) Two bubbles, full linearized gravity approximation 

For two bubbles in the full linearized-gravity approximation, we can derive formulas 
analogous to those used with the scalar field evolution. As in the quadrupole case, let the 
bubbles be nucleated at t = and at 2; = ±d/2, with r/d = 1.2 and the same time cutoff 
function. The problem possesses rotational symmetry about the 2;-axis, so without loss of 
generality we take ky = 0, = cjsin^, kg = a;cos^. Using the same conventions as in the 
quadrupole case, the stress-energy tensor components are given by 

noo r — a /»27r 



— TT 



(24) 



Ja Jo 

The (f) integral can be done explicitly using the identity 

cos X ^ 2^77 J„(/3), (25) 

resulting in the following expressions: 

poo pTT — Ot 

T^^ (k, aj) = q^ dt e''^H^C{t) / dO sin^ d cos{kgt cos 6 + ^k^d) 
3 Jo Jo 

X [Jo{kxtsm9) — J2{kxtsm6)] , (26a) 

Tyy{k, uj)^^ dt e'^^H^Cit) / de sin^ cos{k^t cos 6 + ^k^d) 
3 Jo Jo 

X [Jo {ka;t sin e) + J2 {ka;t sin 9)] , (266) 

T,,(k,u;) = dte'^'H^Cit) ddsin9cos^9 
3 Jo Jo 

X cos(/cat cos 6' + ^kzd)Jo{kxt sin 9), (26c) 

Txz{k,Lo) = — ^ dte'^^H^Cit) d9sin^9cos9 
3 Jo Jo 

X sin(A;^tcos^ + ^kgd)Ji{kxt sin 9). {2Qd) 

Note that Tj;y = Tyz = 0. The energy radiated in gravity waves simplifies to 
dE 

= Gu^ T,,{k, u) sin^ e + T^^(k, u) cos^ ^ 

2 

— Tyy(k, w) — 2Ta;2(k, ct)) sin^cos^ . (27) 
9 



In Fig. 4 we compare the envelope approximation with the previous results using scalar- 
field evolution. As in the quadrupole case, the agreement is excellent, with the envelope 
approximation power being slightly greater. 



III. Numerical Methods and Results 



(a) Bubble nucleation [7] 



The envelope approximation closely reproduces the gravitational radiation from a two- 
bubble collision, even showing the same features in the spectrum. This gives us confidence 
to apply the approximation to the situation where many bubbles are nucleated, collide, 
and transform all of space to the true vacuum. In a phase transition the bubble nucleation 
rate per unit volume, F, is in general a function of time, due to its dependence upon the 
temperature of the universe or the evolution of other fields. Since very generally F is the 
exponential of some action, we write it as 



The prefactor C is expected to be of the order of M.^., where M is the mass or energy 
scale characterizing the transition. The tunnelling action A{t) must be greater than order 
unity; otherwise, the transition will not proceed via bubble nucleation but through spinodal 
decomposition, because of the very small potential barrier between the false and true 
vacuum states. (For the form of the potential in Eq. (5), A 3> 1 obtains for e -C 1.) The 
fact that the nucleation rate varies with time is crucial to the completion of the phase 
transition; moreover, how fast it varies with time determines the distribution of bubble 
sizes. As a rough rule, the phase transition completes when one bubble is nucleated per 
Hubble volume per Hubble time, i.e. when T{t)/H^ ~ 1. Denote the completion time, 
about which we shall be more specific, by t*- Expanding A{t) about gives 



where A^ = A{t^). In any sensible model, /3 > 0; ie., the nucleation rate grows with time. 
As we shall see, (3~^ sets the time/length scale for the phase transition. 

We now derive some important results for the exponential nucleation rate. The fun- 
damental quantity is p{t), the probability that a given point in space remains in the false 
vacuum at time t. It is given by p{t) = exp[— /(t)], where I{t) is the expected fraction of 
space occupied by true-vacuum bubbles at time t, without regard to bubble overlap: 



where a{t) is the cosmic scale factor, r(t, t') is the coordinate radius at time t of a bubble 
nucleated at time t', and to is the time at which the phase transition begins [12]. For 



r{t) = Ce-^(*) . 



(28) 





(29) 




(30a) 
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simplicity and consistency with our previous neglect of the expansion of the universe, 
we take a to be constant and or(t, t') — t — t' . The neglect of the expansion is justified 
provided the duration of the transition is less than H~^. The second assumption implies 
that bubbles expand to a size far greater than that when nucleated, which is well-justified 
in the cases of interest. Taking ^ —oo with little error, it follows that 

m ^ ^r(t). (306) 

Then the false-vacuum fraction is given by 

p{t) = e-^(*) ^ e-8-rW/^^ (31) 

where the exponentiation accounts for the bubble overlap. From p{t) we can compute the 
duration of the phase transition and distribution of bubble sizes. The start and end of the 
transition are somewhat difficult to define precisely, but this ambiguity is not important. 
To be specific, we can define the start of the transition to be the time tm when p{tm) = 
e~'^ ^ 1, i.e.,m <^ 1. Similarly, we define the end of the transition to be the time when 
p(t*) = 0, i.e. M » 1. The duration of the transition is thus 

5t = U-tm = ln(—] (3-\ (32) 



and depends only logarithmically upon the precise definition of the start and end of the 
phase transition. Note that the duration of the phase transition is set by and that 
for consistency our neglect of the expansion of the universe requires the duration to be less 
than a Hubble time: (3~^ <C H~^. 

The density (per unit volume) of bubbles of a given radius r at time t is related to T{t) 
and p{t) by 



t'=t-r 



~ ^^exp[-/(t)e-^'" - /3r]. (33) 
The distribution of bubble sizes attains its maximum at 

r-(t) = ^ln/(t) (34) 

and has a width of order f3~^. In discussing gravitational wave production it is more 
appropriate to examine the energy-weighted bubble distribution. Since the energy carried 
in the expanding wall of a bubble is proportional to its volume, this distribution is obtained 
by multiplying dn/dr by r^: 

{^!^)=?^^r^^[-me-^'-Pr]. (35) 

11 



This distribution is peaked at a radius twice as large as dn/dr (see Fig. 5). 

Finally, how is the key parameter /5 related to H and Ail Since sets the scale 
for all time evolution in the universe, on very general groimds we expect /3 = —{dA/dt)t, 
to be of the order of A{t^) / H~'^ , or ~ H~^/A^. If the transition is to proceed 
via vacuum bubbles, A* must be much greater than one, so the assumption that the 
transition is "fast", /3~^ <^ H~^, should generally be satisfied. We can also estimate A*: 
r(t*) = M'^e~^* ~ ~ Al^/mpi*, which implies that should be of order ln(mpi/A^). 

(h) Numerical results 

In the previous sub-section we have motivated the use of an exponential nucleation 
rate; specifically, 

m = Toe^K (36) 

As we have discussed, sets the fundamental time/length scale: both the duration of 
the transition and the typical bubble size are of order a few We use a spherical 

volume, and choose Fq so that on average, the desired number of bubbles are nucleated in 
the sample volume. Our main calculations (see below) use a sample volume with a radius 
of 4.46/?-^ and Fq = 1.38 x IQ-^P^. These parameters yield an average of around thirty 
bubbles in the sample volume. For our five different nucleation runs, the time at which the 
phase transition completes varied from to with an average of 5.63/?"^. (We 

use the time when the last bubble is nucleated in the sample volume as the completion 
time of the phase transition.) 

A "nucleation run" proceeds as follows. For each time step we: (1) Perform a Monte 
Carlo integration over the sample volume to determine the fraction still in the false- vacuum 
state. This takes into account all true-vacuum bubbles nucleated prior to the time step. 
(2) Multiply T{t) by the sample volume still in the false vacuum and length of the time 
step to get the mean number of bubbles nucleated during the time step. (3) Use this mean 
number to generate a Poisson-distributed random number of bubbles nucleated during the 
time step. (4) Nucleate this number of bubbles at random points in the false-vacuum and 
during the time step. We choose the time step so that the time from the first nucleation 
to the disappearance of all the false vacimm is around 100 time steps. 

The initial conditions for the simulations consist of nucleation sites and times within 
the spherical volume produced by the above procedure. We use a spherical boundary 
to minimize edge effects. We also utilize reflecting boundaries: When a bubble expands 
into the boundary, the part of the bubble that passes the boundary is ignored. In the 
context of the envelope approximation, this is equivalent to having another bubble outside 
the boundary expand into the bubble in question, with their first contact occurring at 
the boundary surface. Note that in the envelope approximation, refiecting and absorbing 
boundaries are equivalent. 

The aim of the simulation is to numerically evaluate Eq. (14) for the stress-tensor 
components Ty(k, a;) for the given initial conditions, and then Eq. (15) for the spectrum 
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of gravitational radiation produced. Note no radiation is produced before the first bubble 
collision occurs or after the transition to true vacuum has completed; thus these are the 
limits for the time integration. Our numerical calculations proceed as follows: First we 
choose values for uj and the direction k. Then we numerically evaluate the integrals 
in Eq. (14). The time partition is chosen for around 100 time steps, and the angular 
partition is chosen to give 20 divisions per wavelength at a given radiation frequency. These 
partitions produce numerical integration results accurate to within 5% at all frequencies 
of interest. 

Since bubble nucleation is an inherently random process, we want to find the power 
spectrum of gravity waves, averaged over many different nucleation realizations. We have 
calculated the radiation from five simulations averaging 30 bubbles each, ranging from 17 
to 38 bubbles. Figure 6 shows a cross-sectional slice through the equator of the spherical 
volume for one nucleation run; the combined bubble envelope, which is the radiation 
source, is shown at several times. For each simulation, we have computed the radiation 
in six directions, along the ±x, ±y, and ±z axes. The results for dE/dwdVL for several 
representative directions are displayed in Fig. 7. 

Several features are evident from these results. First, all of the spectra peak at a 
characteristic frequency which roughly corresponds to the size of the largest bubbles at 
the end of the transition. Second, the power generally increases as the number of bubbles 
in a given volume decreases, as we expect {dE/dcvdQ should vary as 

Ar-5/3). Finally and 

most obviously, the various runs and directions of observations have very large fluctuations 
in both total power and spectrum shape. This is because with the relatively small number 
of bubbles, the source has large inhomogeneities on the scale of the sample volume. One 
direction in one simulation {i.e., one out of 30 probes of dE/dwdVt) exhibits a "beaming" 
effect, with ten to a hundred times more power than average at frequencies higher than 
several times the peak frequency. 

Figure 8 shows the power radiated in gravity waves per octave averaged over all five 
simulations, six directions per simulation (excluding the one anomolously "hot" direction); 
the error bars indicate the statistical deviation of the mean. On the same plot we show our 
previous results for the collision of two scalar- field bubbles [6] . Recall the previous results 
depend upon the cutoff time r, which corresponds to the duration of the transition. To 
compare the old results with the present ones, we set r — 3.2//3. We choose this value for 
T because it is the mean energy- weighted bubble size at the end of the phase transition (c/. 
Fig. 5). We have divided our results by the total vacuum energy of the sample volume, so 
the results shown are the fraction of false-vacuum energy liberated in gravitational waves 
per octave. The two-bubble and many-bubble results are remarkably similar: behavior 
at high and low frequencies is almost identical, and the peaks are at almost the same 
frequency. The overall normalization of the many-bubble case is higher, by a factor of five 
or so on the low-frequency side and by an order of magnitude on the high-frequency side. 
This increase is not unexpected: In the many-bubble case, each bubble collides with many 
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others, increasing its total radiation. The excess of high frequency power also makes sense, 
since the size of a given bubble roughly determines the frequency at which it radiates. In 
the two-bubble case, both bubbles are the same size; in the many-bubble case, smaller 
bubbles are nucleated late and increase the high-frequency power. 

The total fraction of vacuum energy released in gravitational radiation is computed by 
integrating dE/dw and dividing by the total vacuum energy released: 



The peak of the energy spectrum is given by ct^max ~ 1-6/9. We note that the charac- 
teristic frequency defined by the maximum of the energy-weighted bubble distribution 
a;* = I'K jfE ~ 1-6/3 coincides with o^max- Comparing these results with our previous two- 
bubble results (taking r = 3-2//3), the peak of the spectrum occurs at the same frequency, 
while the total fraction of energy liberated in gravity waves is about a factor of five larger. 

To verify that our results are not dominated by edge effects, we ran one large calculation 
in a spherical sample volume with twice the radius of the above simulation and the same 
nuclcation rate, 180 bubbles in all. The radiation spectrum dE/duj, integrated over the 
six observation directions, is plotted in Fig. 8. The shape of the spectrum is close to the 
average of the previous cases, and Eow/Eyg^c — 0.47Gpvac//5^- The close correspondence 
(only 5% difference) with our smaller simulations demonstrates that edge effects are not 
important. 

IV. Statistical Bubble Approximation 

Now we discuss two approximations for computing the gravitational- wave production 
from the collision of many bubbles as the incoherent sum of the radiation from individual 
bubbles. To this end, we use the multipole-radiation formalism and the envelope approxi- 
mation applied to a typical bubble of size R, and integrate over the distribution of bubble 
sizes. The advantage of such an approach is ease of computation: We can immediately 
calculate the radiation spectrum for a given nucleation rate, without recourse to a many- 
bubble simulation. Of course, this approach neglects coherent effects between bubbles and 
is less accurate than our calculations of the previous Section, but as we shall see it still 
gives a reasonable estimate of the gravitational radiation produced. 

Consider a single expanding bubble in a sea of bubbles. In the envelope approximation, 
when the bubble in question begins interacting with other bubbles, portions of the bubble 
surface are "eaten." Neglecting the interaction regions between bubbles, the stress tensor 
for what remains of the bubble is 




(37) 




(38) 



where is defined as 




1, surface in direction Q remains uncoUided 
0, otherwise. 



(39) 
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When the bubble wall has completely disappeared due to collisions with other bubbles, the 
function Q is zero. This means that an individual bubble can be treated as a source that is 
bounded in space and time: At, Ax < 0{t), where r is the duration of the phase transition. 
The size and energy density of a given bubble at any time during a phase transition is 
known; the function varies from bubble to bubble. Our statistical approximation boils 
down to estimating an "average" 0(0, t). 

(a) Multipole-radiation formalism 

Using the tensor spherical harmonics presented in the Appendix, we can expand the 
transverse-traceless part of the metric perturbation {i.e., the gravity wave piece) in the far 
field zone (r 3> r) as follows [13]: 

I 



1=2 m=-l 



^/-^{t - r)Tf ''-(O) + |^5^-(t - r)T- 



(40) 



where the /'"^ and the are the "mass" and "current" multipole moments of the source, 
respectively, defined by 

^/'^(t) = 8(-z)'+2 j r'^dr'dn'dt'dw e-'^^'-''^Tpg{t' , r', n') 

X a_2(0Tjr'''"^(^^')*J/-2(a;r') -ao(OTj^'''-(0')*jKa;r') 



(41a) 



dt^ 



.-(.) = 8(-0-/.-™^e-M.-O.,^(,,.^^ 



(416) 



The tensors T^*''"^(n) and the coefficients are given by Equations (A. 2) and (A. 4) in 
the Appendix, ji is the spherical Bessel function of order I, p and q are indices of the 
tensor components, and Tpg is the sum of the stress-energy tensor of matter [Tpq) and the 
Landau-Lifshitz pseudotensor for the effective stress-energy of the gravitational field. Since 
we are using the linearized gravity (appropriate since all gravitational effects are weak; see 



Ref. [6]), Tpq 



Tpq. The effective stress-energy tensor for gravity waves is 



327rG' 



TT\ 
ij I 



(42) 



and the power radiated in terms of the multipoles: 

G 



l.m 



Jl+1 

" Tim 


2 


" aim 


+ 


dt^+^ 




dt^+^ 



(43) 
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where ( • • • ) indicates a spatial average over several wavelengths. 

The expressions (41) can be simplified considerably. First, only the exponential factor 
and the Bessel functions depend on a;. The lo integral can be performed explicitly using 
the identities 



<i; 

>1, 



(44) 



where Pn are the Legendre polynomials. This gives 



ji foo r r 

|^/'-(t) = 87r(-l)'y^ r''dr'J dr/Tp,(t - V, /, J]') [«-2(0Tj^''^"^(^^')*^'/-2(r7) 



pq 



(45) 



Now for a single bubble nucleated at t = 0, the 0(3, 1) symmetry of the scalar- field bubble 
configuration allows transformation of the radial integral: 

i=-2,Q,2 \ V J 



where 



A'''^-(n') = a,,_K/)T5['^-(n')*ip(0')x,(0'). 



, (46) 



(47) 



Here s = a/t*^ — (t — rrj)"^ is the only quantity the scalar field depends upon, 0'(s) = 
d(p/dr{r,t = 0), and x is the unit vector in the direction of Q'. Explicit expressions for 
the A'' ''"^ are given in the Appendix. The further substitution 



t-y 



gives 



ds f dip 



ds 



i=-2,0,2 



t/2 y + s 



t-y 



dn'A^+'^^'^{n')Q{Q',y) 



(48) 



This expression can be simplified by noting that (1) d4>/ds is non-zero only when s is 
smaller than the initial bubble radius, and (2) the angular integral is zero until the bubble 
first collides, which by assumption is only after the bubble has expanded by a very large 
factor. Thus y ^ s and to a very good approximation s can be set to zero in the integrand 
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of the y-integral. Now using Eq. (9b), the s-integral can be performed to give Pvac/3; 
Eq. (48) becomes 




it) = ^(-l)Vvact^ / 




(49) 



where 



(50) 



Note vanishes (by spherical symmetry) until the bubble first collides; likewise, 

Qi+i,im vanishes for t > where R is the size of the bubble when its surface has completely 
collided, since 0(0, t > -R) = 0. Analagous formulas holds for the "current" multipole 
moments S''"*, but they vanish because the tensor-spherical harmonics contracted with the 
unit vectors are identically zero (see Appendix). 

The multipole radiation from a single bubble is only nonzero for < t < 2R. This 
makes sense physically. The bubble wall propogates outwards at essentially the speed of 
light; even though the bubble does not begin to radiate until it first collides, the first 
radiation from the bubble still reaches an observer at distance r at time t = r. Likewise, 
radiation generated on the opposite side of the bubble from the observer will arrive at time 
t — r + 2R, since the diameter of the bubble when it disappears is 2R. 

Equation (49) provides the key to computing the gravitational radiation from a single 
bubble. The evaluation of SP'^ /dt^ only involves computing 6'"''*'''"(t), which depends 
upon the "collision history" of a given bubble. We present two different estimates for 
Ql+i,lm{ty The first is analytical, based upon the fraction of the bubble surface that 
remains uncoUided at time t; the second estimate is derived from our numerical simulations 
of bubble collisions. 

Once we have d^P'^ /dt^ in hand for a single bubble, it is a simple matter to calculate 
the power and the spectrum of the gravitational radiation by summing incoherently over 
the distribution of bubbles: 



where ('^/'"^(a;, i?) is the Fourier transform of d^P'^ /dt^{t, R). The R in the argument 
refers to the multipole radiation from a single bubble whose size at the end of the phase 
transition is i?, and dn/dR is the distribution of bubble sizes, cf. Eq. (33). Note that the 
expressions in Eqs. (51) and (52) are for power and energy per unit volume (since dn/dR 
is the bubble size distribution per unit volume). 




(51) 



(52) 
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(b) Analytic approximation to 0'+*'''^(t) 

First we consider a very simple analytic approximation. Recall that 0'+«'''^ is an 
integral over the uncoUided bubble envelope. It must depend upon the fraction of the 
bubble wall that remains uncoUided at time t. This fraction, f{t, tn), is given by [7] 

/(t,tfl) = e-^«+^(*«) (53) 

where tu is the nucleation time for a bubble which has radius R at the end of the phase 
transition, i.e., tR = — R. Since 0'"'"*'^'"(t) must vanish at early times when the bubble 
has not yet collided (/ = 1), and at late times when the bubble surface is completely 
collided (/ = 0), Q^+*'^"^(t) can be expressed as a sum of the terms (1 — /)"^/"' (n, m = 
1, 2,. . . ). We take as a simple ansatz 

e^+M-(t) = c/(l-/) (54) 

where c is an undetermined normalization constant. 

We now estimate gravitational radiation using the quadrupole {I = 2) term of the 
multipole expansion. Taking only the simplest {i = —2) term of Eq. (49), we approximate 
/2'"(t) as 

|^/^-(t, R) = y Cpvac«' J f{u, R) [1 - f{u', R)] du' (55) 

where we set = as the time origin so that Ir = — i?, u = t — tii = t + R, and T oc e^*, 
giving 

f{u, R) = exp [-Me^(^-^^ + Me"^^] ; (56) 

recall M = /(t*). Note by setting = instead of tn = 0, and by using u = t — tn, 
we take into account that bubbles of different sizes are nucleated at different times. The 
power radiated in gravitational waves per unit volume is given by the incoherent sum over 
the distribution of bubble sizes, Eq. (51); we have replaced the sum over m in Eq. (51) 
with a factor of 5. Likewise, the energy spectrum of gravitational waves per unit volume 
is given by Eq. (52). 

In Fig. 9 we show P{t) and in Fig. 10 we show udEay^ /du. At low frequencies, the 
spectrum behaves just as in our many- bubble simulation; however, at high frequencies 
it falls off more rapidly, and the peak of the spectrum is about a factor of three lower. 
The deficiency in high frequency power traces to the fact that we have neglected the 
sharp cusps which form in bubble collisions. In this approximation the fraction of vacuum 
energy released in gravity waves is Eqw/Evac = 1.5c^{H/P)'^, so that c ~ 0.2 reproduces 
the result of our many-bubble simulations. We have tried a range values for powers of 
/ and 1 — / in Eq. (54); these alternatives change only the overall normalization of the 
spectrum (decreasing with increasing powers). 
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(c) Numerical estimation of 



It is straightforward to extract from our numerical simulations average multipole mo- 
ments for a single bubble, and thereby determine a normalized approximate spectrum. 
Specifically, for 70 individual bubbles, nucleated in two of our five smaller simulations, we 
have computed d} {t) / dt^ by using Eq. (49), for Z = 2 and I = 3 (the quadrupole and 
octupole moments). These multipole moments give the gravitational waveforms from a 
given bubble via Eq. (40); quadrupole and octupole waveforms for a representative bubble 
are displayed in Fig. 11. The radiated power P{t) for the same bubble is shown in Fig. 12; 
note that the octupole power is only around 10% of the quadrupole power, so we can safely 
assume that contributions from higher moments are negligible. 

To calculate an "average" energy spectrum from the 70 bubbles, we first construct 
scaled multipole moments, removing a factor of R from the time variable and a factor of 
Pvac-R^ from d^P^{t)/dt^ for each bubble. Then we calculate the energy spectrum for each 
bubble using these scaled moments, and average over the 70 bubbles to give the average 
spectrum for a bubble. This average energy spectrum, shown in Fig. 13, is quantitatively 
similar to the two-bubble spectrum, peaking at around ujR = 4.6 with roughly the same 
overall normalization, but dropping off much faster on the high-frequency side of the peak. 
We then follow the same procedure as in the previous subsection, integrating over the 
bubble-size distribution to give the energy spectrum per unit volume of the radiation from 
the phase transition. The energy per octave is compared with the spectrum from our 
many-bubble simulation in Fig. 10. 

Overall, this statistical approximation does a reasonable job. As with the previous 
analytical approximations, it closely reproduces the low-frequency behavior. The peak of 
the approximate spectrum has about the correct amplitude, though the peak frequency 
is low by about a factor of two. The most obvious discrepancy is again the rapid high- 
frequency drop of the approximate spectrum, which falls off as for large frequencies, 
in marked contrast to the many-bubble calculation, which falls off as As before, the 

deficiency of high-frequency power is due to the neglect of cusps. In this approximation the 
fraction of energy radiated in gravitational waves is 

Eow/E^ac = 0.036(if//3)2, compared 

with 0.06(if//3)2 in the many-bubble simulations. 

The utility of our pair of statistical approximations lies in their computational ease. 
The difference in computing time between these approximations and the many-bubble 
simulations is enormous. In Fig. 10, the analytical statistical approximation required 
negligible computing time, the numerical statistical approximation used around an hour of 
computing time, while the many-bubble points required several weeks on the same machine. 
Furthermore, our averages for the radiation from a single bubble can be used with any form 
for the bubble nucleation rate, by substituting the appropriate form for dn/dR in Eq. (52). 
These approximations give a rough, but quick, estimate of the gravitational radiation from 
any strongly first-order phase transition for which the bubble nucleation rate is known. 
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V. Discussion and Concluding Remarks 

Before summarizing the present work, let us place it in context by reviewing our pre- 
vious work. Based largely on dimensional estimates, it was argued that the gravitational 
radiation produced by the collision of vacuum bubbles in a strongly first-order phase tran- 
sition could account for a substantial fraction of the vacuum energy released [5]. This 
conjecture was verified in our previous numerical work [6] where we calculated the grav- 
itational radiation resulting from the collision of two vacuum bubbles by evolving the 
scalar-field configuration corresponding to two vacuum bubbles nucleated simultaneously 
and separated by distance d. This calculation was carried out in the linearized-gravity 
approximation and the expansion of the universe was neglected, both assumptions being 
valid provided that the duration of the transition is less than a Hubble time. We found 
that the amount of radiation emitted is indeed significant and only depends upon the dura- 
tion of the collision — and not the fine-scale details of the bubbles. The fraction of vacuum 
energy liberated into gravitational waves is ^^owZ-^'vac — 1-3 x 10~^ {HtY, valid for t 
where r is the duration of the transition. Unfortunately, this work depended upon the 
phenomenological parameter r; moreover, it is a bold extrapolation to use the collision 
of two bubbles to model a realistic phase transition, which consists of many bubbles of 
different sizes colliding. As noted earlier, it is beyond present computational capabilities 
to collide more than a few bubbles by scalar-field evolution. 

These drawbacks led to the present work: the development of a workable approxima- 
tion to study the gravitational radiation from hundreds of colliding bubbles. Motivated 
by the fact that our two-bubble results only depend upon the gross features of the colli- 
sion, we developed the envelope approximation described in this paper. In the envelope 
approximation an expanding bubble is treated as a very thin shell of energy (equal to the 
vacuum energy it liberates); when bubbles collide only their envelope is followed and their 
overlap (interaction) regions are ignored. By considering the collision of two bubbles we 
showed that the envelope approximation accurately reproduces our previous two-bubble 
results; e.g.^ the energy spectrum in gravitational waves agrees to around 20%. 

Having established the validity of the envelope approximation, we nucleated hundreds 
of vacuum bubbles in spherical volumes with a nucleation rate that grows as e^^: specifi- 
cally, 127 bubbles total in five small simulations and 180 bubbles in one large simulation. 
(It is argued in Ref. [7] that such a functional dependence for the bubble-nucleation rate 
applies with great generality.) Using the envelope approximation we computed the frac- 
tion of vacuum energy released in gravitational waves and found: Egy^/E^q_c — 0.06 (i?//?)^ 
with the energy spectrum peaking at a frequency uj ~ 1.6/5 {H'^ = 87rGpvac/3). With this 
nucleation rate the duration of the phase transition r ~ 3//?; using this fact, it follows 
that the fraction of energy liberated in gravitational waves is about five times the estimate 
based upon our previous two-bubble results, with the spectrum peaking at about the same 
frequency. We believe that additional energy is released in gravity waves because a given 
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bubble collides with many other bubbles rather than a single bubble. 

In the present work we also have developed two statistical approximations that allow 
simple analytical or semi-analytical approximations to the energy spectrum radiated in 
gravitational waves. Both approximations provide better than order-of-magnitude accu- 
racy and greater ease of calculation, and are particularly well suited to computing the 
gravitational radiation for an arbitrary bubble nucleation rate. 

On very general grounds it has been argued that which controls the time/length 
scale of the phase transition, is of the order of a few percent of or greater (see Section 
IV or Ref. [7]), indicating that the fraction of vacuum energy liberated in gravitational 
waves in a phase transition that proceeds through the nucleation and collision of vacuum 
bubbles is of order lO""^ or so. We have addressed the potential observational consequences 
of our results in a Letter [8] ; very briefly, in terms of the temperature of the universe after 
the phase transition, the fraction of critical density contributed by gravitational waves 
produced is Oqw ~ 10~^ with characteristic frequency / ~ 10~^ (T/ GeV) Hz. There we 
also discuss the prospects for the detection of such a stochastic background of gravitational 
waves with the coming generation of laser interferometer gravity- wave observatories [14] . 

Two key assumptions underlay all of our work: (1) the use of linearized gravity and the 
neglect of the expansion of the universe; and (2) the assumption that the bubbles expand at 
constant acceleration (put another way, all the vacuum energy liberated is converted into 
the kinetic energy of the bubble wall). As we discussed in Ref. [6] the first assumption is 
valid so long as the duration of the phase transition r ~ is much less than the Hubble 
time H~^. This should be satisfied for most phase transitions as is expected to be 
only a few percent of H~^. However, there are situations where this condition may not 
be satisfied, e.g., in some models of extended infiation [15]; moreover, such situations are 
very interesting since our results indicate that the fraction of energy liberated in gravity 
waves approaches unity. We are currently trying to generalize our results by relaxing the 
first assumption [16]. 

The second assumption is that all the vacuum energy liberated by a bubble goes into the 
kinetic energy of its wall. This is certainly true for a bubble nucleated at zero temperature, 
but it may not be a good approximation to one nucleated at finite temperature because of 
the interaction of the bubble with the ambient thermal plasma. The second assumption is 
certainly well justified for models of first-order infiation where the universe has undergone 
extreme supercooling during the infiationary epoch so that the temperature of the Universe 
when the phase transition occurs is exponentially small. Whether or not this assumption 
applies in a first-order phase transition that only undergoes moderate supercooling is an 
open question. In this case it is not implausible that much or even most of the latent heat 
released is dissipated into heat rather than the bulk motion of the expanding bubble front 
(here wc have used the term latent heat rather than vacuum energy). The motion of a 
bubble wall in this circumstance is not a simple matter to analyze: both the microscopic 
interaction of the ambient medium with the bubble front and bulk hydrodynamics are 
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important. Though much work has been done the results are not conclusive [17]. 

The strength of a first-order phase transition can be characterized by the ratio of 
the latent heat (per unit volume) released to the energy density of the ambient plasma, 
given by the fourth power of the temperature at which bubble nucleation commences: 
7 = Pvac/T^nuci ^otc that the increase in entropy per comoving is proportional to 7^/^. For 
an inflationary transition 7 ^ 00, while for a weakly first-order transition 7 is of order unity 
or less. For very large 7 it seems clear that the bubbles must behave as vacuum bubbles 
(all the latent heat liberated goes into accelerating the bubble wall). What happens for 
moderate values of 7 is still unclear and probably depends on the specific phase transition 
under consideration. The latent heat could simply be dissipated viscously, in which case 
little gravitational radiation would be produced; or the latent heat could be converted into 
the bulk motion of the fiuid (at some velocity less than the speed of light) , in which case 
an appreciable amount of gravitational radiation could still be produced. This important 
issue is currently under study [18]. 
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Appendix: Tensor Spherical Harmonics 

In this appendix we review the formalism of tensor-spherical harmonics, and calculate 
the tensor contractions needed in Sec. IV. Generally, we follow the notation presented 
in [13]. A set of basis vectors ^"^ are coupled to form the traceless and symmetric basis 
tensors t™: 

+1 +1 

t"^= ^ J2 (llm'm"|2m)e"''®^'"" (Ala) 

m' = — l m"= — l 

where {h TO2K3TO3) is the Clebsch-Gordan coefficient for adding angular momenta h 
and I2 to obtain ^3. In terms of the Cartesian basis vectors Bx, ey, and e^, these symmetric 
basis tensors are 

1 i 

t^^ = -{Bo: ^e^-Gy ey) ± -{e^ By + By B^)] (A.I6) 

1 i 
t^^ = T-i^x ® + e^ ® e^) - -(e^ O e^ + e^ ® e^); (A.lc) 

t° = ^{2Bz <S> Bz - Bx <S> Bx - By <S> By). (A.ld) 
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Then the relevant tensor spherical harmonics are 



I' 2 

rp2r,/m^ J2 J2 a'2m'm"|Zm)y''^'t^", (A2) 



m'=-l' m"=-2 



where I' = I, Z± 1, Z±2. This represents the combination of an orbital angular momentum 
I' and a spin angular momentum 2 to give total angular momentum I. These spherical 
harmonics are eigenfunctions of the orbital angular momentum operator with eigenvalue 
l{l + 1), like the more familiar y'"^. Also, we have the "pure spin tensor harmonics" for 
spin 2, 

rj.E2,im ^ ^^^i^rj.2i+2,im ^ ao{l)T^^'^'^ + a_2 (OT^ (A3a) 

where 



"^(^)"^/ 2(2/ + 1X2] + 3) ' (^-^"^ 



= (AAb) 



-°(^)"^/ (2/-i)(2/+3r ^^-'^^ 



/,^ / (/ + l)(^ + 2) , , , , 

-^-^^^^"^/ 2(2/ + l)(2/-l) - ^^-^^^ 

Under rotations around the radial vector, the harmonics (A. 3) transform like the compo- 
nents of the polarization tensor of a pure spin-2 state, but they are not orbital angular 
momentum eigenfunctions. 

Using the above definitions, it is straightforward to calculate the explicit forms for the 
necessary tensor contractions. Equation (47) defines the auxiliary quantity 

where x is a unit vector in the direction fl. The expressions for m < follow from the 
identity 

r£2l',lm ^ -^y' +l+m ^1-^2 I' ,1 —m^* 

For I = 2 (quadrupole moments): 
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^0,22^ 1 /Ae-2i^sin2^ (A7) 

4 V 57r ^ ^ 

^2,22 ^ _ 5 /|^^_2^<^^^ ^ ^.^2 ^ ^ 3 sin^ ^) (A.8) 
28 V OTT 

^4,22^1 /^g-2^^g.^2^ 

28 V Stt ^ ^ 



^0,21 = _i J Ae-»<^ sin 26 (A.IO) 
4 V 57r ^ ^ 

l^'^^ = ^y^e-*'^(8-2cos2^-sin4^) (^.11) 



1 / 3 




^*'''° = ^\/^(l+3'=os2e) (A15) 

28 V OTT 

For Z = 3 (octupole moments): 

^/',3m^Q £qj. ^/^2, /' = 4 (A16) 

^1,33 ^ _1 /Ae-3#sin3^ M.17) 

^3,33 ^ 1 /Ae-3#sin3^ (A18) 

18 V Ttt ^ ^ 

^5,33 ^ _1 /Ae-3^<^sin3^ M.19) 
36 V Ttt 

4 V Ttt ^ ^ 

^3,32 ^ _Z /^e-2i<^cos^sin2^ (A.21) 
6 V Ttt 

^5,32 ^ 1 . / Ae-2»<^ COS ^ sin^ ^ ( A.22) 

12 V Ttt ^ ^ 
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^i.3i = -,/-^e-*'^sin^ 



^3,31 




7- 



21 7 7 

^ cos 2^ - — cos 4$ 

80 12 96 




77r 



^5,31 



^1,30 ^ 
^3,30 ^ _ 
^5,30 



> i n I J 



1 5 

- H cos 29 

8 24 



2l7r 



2l7r 



2l7r 



cos 9 



3 3 
- cos 29 

4 20 



^cos^+^cos30 

1 5 

- cos 9 + — cos 3^ 

8 24 



(yl.23) 
(yl.24) 

{A.25) 

{A.26) 
{A.27) 
{A.28) 
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Figure Captions 



Figure 1: The kinetic (dashed curve), gradient (dot-dashed curve), and total (sohd 
curve) energy of a single expanding vacuum bubble. The scales are arbitrary; the bubble 
has initial radius of about 10, in the same units as t. Note that hy t = 20, when the 
bubble has approximately doubled in size, the total energy scales almost exactly as t^, the 
vacuum energy liberated by a bubble expanding at the speed of light from zero initial size, 
and resides equally in the kinetic and gradient energies of the bubble wall. 

Figure 2: A schematic picture that illustrates the envelope approximation. The dark 
lines are the bubble walls, expanding at the speed of light. The shaded areas are the 
interaction regions; the envelope approximation neglects the interaction regions and takes 
into account only the bubble envelopes. Snapshot (b) is at a somewhat later time, and 
three new bubbles have been nucleated. 

Figure 3: The energy spectrum of radiation from two colliding bubbles, in the 
quadrupole approximation. The units are the same for both curves, but arbitrary. The 
dashed curve is the result from detailed scalar field evolution (Ref . [6] ) , and the solid curve 
from the envelope approximation. 

Figure 4: The energy spectrum from two colliding bubbles in the full linearized- 
gravity approximation. The units are the same as in Fig. 3. The solid line is the envelope 
approximation, which reproduces well the results from detailed scalar-field evolution [6], 
the dashed curve. 

Figure 5: The distribution of bubble sizes, both unweighted and weighted by the 
bubble's total energy, for the exponential nucleation rate T{t) = Foe^*. The energy- 
weighted distribution peaks at a radius that is almost twice as large. 

Figure 6: A slice through a spherical sample volume at different times. The volume 
boundary is the outer circle; the bubble walls are the darker curves within the boundary. 
A total of 33 bubbles were nucleated in this volume during the phase transition. 

Figure 7: The differential energy spectrum dE/dwdCl, for three orthogonal direc- 
tions in a simulation with 33 colliding bubbles. Note the large variations in the different 
directions. 

Figure 8: The energy spectrum per frequency octave, divided by the vacuum energy 
of the sample volume. The points with error flags are averaged over the five different 
simulations and integrated over six directions per simulation. The error bars reflect the 
standard estimate for the deviation of the mean. The dashed line is the spectrum for two 
bubbles, calculated using scalar-field evolution [6], with r = 3.2/3. The solid triangles are 
the results of the 180-bubble simulation (integrated over six directions). 

Figure 9: The radiated power P(t) per unit volume in the first statistical approxima- 
tion. The phase transition completes ai t — 0. 
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Figure 10: The energy spectrum as computed in the first statistical approximation, 
with c = 1 (sohd fine) and the spectrum from the second statistical approximation (dashed 
line), compared with the results of our many- bubble simulations (points with error flags). 

Figure 11: The quadrupole and octupole moments for a given bubble nucleated at 
time zero. R is the final bubble radius, (a) The real and imaginary parts of the quadrupole 
moments for m = 0, 1, 2 (the imaginary part of the m = moments vanish), (b) The 
same for the octupole moments for m — 0, 1, 2, 3. 

Figure 12: The quadrupole (upper curve) and octupole (lower curve) contributions 
to the total power radiated from the bubble in Fig. 10. 

Figure 13: The mean energy spectrum of a "typical" bubble, derived by averaging 
over 70 individual bubbles. The standard estimate for the deviation of the mean is around 
±10%. 
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